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ABSTRACT 

We investigate the uncertainties affecting the temperature profiles of dense cores of interstellar 
clouds. In regions shielded from external ultraviolet radiation, the problem is reduced to the 
balance between cosmic ray heating, line cooling, and the coupling between gas and dust. We 
show that variations in the gas phase abundances, the grain size distribution, and the velocity 
field can each change the predicted core temperatures by one or two degrees. We emphasize the 
role of non-local radiative transfer effects that often are not taken into account, for example, 
when modelling the core chemistry. These include the radiative coupling between regions of 
different temperature and the enhanced line cooling near the cloud surface. The uncertainty of 
the temperature profiles does not necessarily translate to a significant error in the column density 
derived from observations. However, depletion processes are very temperature sensitive and a 
two degree difference can mean that a given molecule no longer traces the physical conditions in 
the core centre. 



Subject headings: ISM: clouds 
— Radiative transfer 

Introduction 



ISM: molecules — ISM: kinematics and dynamics — dust, extinction 



The gas kinetic temperature T gas is a central 
parameter for the stability of the clouds, their evo- 
lution, and their ch emistry (e.g. iGalli et al.ll2002 ; 
Keto fc Fieldll2005f ). The properties of pre-stellar 
cores are of particular interest because they have 
a direct impact on the ensuing star formation. 
The understanding of the temperature variations 
within the cores is thus crucial for the interpreta- 
tion of observations and f or the development of th e 
theory of star formation (jBergin fe Tafallall2007l) . 

The early models of the cloud cores were spher- 
ically symmetric and isothermal, in the case of 
both stab l e and collapsing cores (e.g. iBonnor 
19561 : IShu 1977 ). The spherical symmetry will 



along the line-of-sight. The dust colour tem- 
perature is known to be a biased estimate of 



the true grain temperatu re Td us t (jShettv et al 



20091; iMalinen et alj|201ll) . The situation is po 



tentially even worse for T gas where, because of 
abundance variations, the observations may rep- 
resent conditions in a small fraction of the source. 
This is particularly true for cold pre-stellar cores 
where many molecules a re severely depleted (e.g . 
Bergin fc Langerl Il997l; iBelloche fc Andrei |2004 



eventually break down because of rotational flat- 
tening or core fragmentation but, in early evolu- 
tionary stages, it is often a reasonable first ap- 



proxi mation (e.g. lEvans et al.l l2005t lAlves et al 
20011 ). The observations of the detailed temper- 
atures structure can be more difficult because 
the measured intensities are non-linear averages 
of the emission originating in different regions 



Whittetl 1201ft iFord fc Shirlevl l201lt IParise et al 
20111 ). There are measurements suggesting that, 
in the centre of a starless core, the ki netic temper- 
ature can decrease dow n to ~6K (jCrapsi et al 
20071 ; lHariu et al] 120081) . However, our knowl- 
edge of the radial temperature variations depends 
heavily on models that also form the basis for the 
interpretation of molecular line data. 

The principles of the thermal balance of dense 
clouds are well understood, the main factors be- 
ing the heating by cosmic rays, the cooling by line 
emission, and the energy exchange between gas 
and dust ((Goldsmith fc Langerl Il978t iGoldsmitb 
2001). Despite the apparent simplicity of the 
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problem, reliable predictions of T gas are not easy 
to obtain because of the uncertainty of the frac- 
tional abundances and dust grain sizes, and the 
potentially complex radiative transfer effects. The 
problem of line cooling has been examined sepa- 
rately, with the Large Velocity Gradient (LVG) 
approximation (jNeufeld. Lepp &: Mehrickl [l ogs') 



and with Monte Carlo methods (Juvelaetal 



20011 ). the latter also enabling the study of the 
role of an inhomogeneous medium. 

In this paper, we investigate the uncertainties 
of the modelled T gas profiles of dense cloud cores. 
The cores are known or are suspected to have sig- 
nificant radial variations in the gas phase abun- 
dances, grain size distributions, and the veloc- 
ity field. We wish estimate the maximum effects 
on T gas that could arise from these systematic 
changes. This is important for the evaluation of 
the uncertainties of chemical model and, more di- 
rectly, the interpretation of any line observations 
of dense cores. 

We describe our models in Sect. [H The re- 
sults are presented in Sect. G2 both for homoge- 
neous models (Sect. 13. ip and for a Bonnor-Ebert 
type cloud (Sect. lX2"1) . and our conclusions are pre- 
sented in Sect. [5j 

2. The Modelling 

We examine the thermal balance of spherical 
clouds without small scale inhomogeneity. All 
models are divided to one hundred concentric 
shells, the innermost shell being 6% and the out- 
ermost shell 0.6% of the outer radius. The same 
discretization is used for in continuum and line 
radiative transfer calculations and for the deter- 
mination of the temperature profiles. To separate 
the radiative transfer effects from those of the den- 
sity, we start by examining homogeneous clouds. 
The dust temperatures are determin ed with Mont e 
Carlo radiative transfer ca lculatio ns ( Juvelal20 05') , 
using the dust model of iDrainel (|2003 ) a nd the 
i nters tellar radiation field (ISRF) given by iBlack 
([19941 ) . The rate for the energy exchange between 
gas and dust is calculated from 

A g . d = 2 x 10- 33 H# 2 )/cm 3 ] 2 (T gas - T dust ) 

{Tgas/lO.OKf^ergcm^s- 1 . (1) 

and the rate for cosmic ray heating from 



(|Goldsmithl l200ll) . Following iGoldsmithl (|2001l ) , 

we calculate the line cooling A gas by 12 CO, 13 CO, 
C 18 0, C, CS, and o-H 2 0, multiplying the last two 
rates by 10 and two, respectively, to take into ac- 
count the cooling by other species. The radiative 
tran sfer problem is solv ed with Monte Carlo meth- 
ods (IJuvela et al.ll200ll). wi th abundances given in 
Table 1 in lGoldsmithl (j200ll ). As part of the Monte 
Carlo simulation, the program saves the net radia- 
tive cooling rates for each model cell. 

As a more realistic cloud model we examine a 
one solar mass, marginally stable Bonnor-Ebert 
sphere with a value of the stability parameter 
£ = 6.5 (jBonnorlll956h . The radial density profile 
is calculated assuming a constant temperature of 
10 K. The variations of kinetic t emperature modify 
the density profile only slightly (|Evans et al 1 120011) 
and this will not affect the conclusions drawn from 
the models. An external UV field can affect the 
cloud temperatures but, in principle, only in a thin 
surface layer. We consider this heating through 
the photoelectric effect in connection with some 
Bonnor-Ebert models. To estimate the photoelec- 
tric heating rate Tpe, we first calculate the energy 
that dust absorbs in the energy range 6-13.6 eV. 
The values are obtained from the cont inuum radia- 
tive transfer calculations. Following IJuvela et al 
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(2003), TpE is obtained by multiplying these en- 
ergies with an constant efficiency of e — 0.029. 
The dust temperatures are calculated assuming 
the clouds are heated by the normal ISRF, even 
when UV heating of gas is ignored. Photoelec- 
tric heating is almost completely eliminated if the 
cloud is surrounded by a dust layer of Ay ~ 2 m 
(see Sect. 13.2. 1|) . This would reduce the central 
dust temperature only by ~0.5 degrees, less than 
the uncertainty associated with the selection of a 
dust model. 

Below we modify several model parameters to 
determine how their variations are reflected on the 
gas temperature. These include (1) the thermal 
coupling between gas and dust, Ag^, that depends 
on the grain size distribution, (2) the abundances 
of the cooling species that depend on the degree of 
depletion, (3) the cosmic ray heating T gaStCr that 
depends on the rate of cosmic rays, and finally (4) 
the large scale infall motion and (5) the small scale 
velocity field that both affect the radiative cooling 
A g . The degree of modification for each parameter 
is indicated in Table Q] and discussed further in the 
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following section. 
3. Results 

3.1. Homogeneous models 

As a first test, before modifying any param- 
eters listed in Table [1] we compared o ur calcu- 
lations with those of lGoldsmithl ( 2001 ') who es- 
timated the line cooling with LVG modelling as- 
suming a velocity gradient of 1 kms -1 pc _1 . The 
cloud is taken to be well shielded from the ex- 
ternal UV field so that the heating through the 
photoelectric effect can be neglected. Correspond- 
ing to the idea of the LVG models, we calculated 
T gas at the centre of homogeneous, microturbu- 
lent spheres with a radius of R = 1 pc and the 
line Doppler width equal to lkms -1 . The corre- 
spondence was found to be good, mostly within 
~10% in T gas (see Fig. QJi), in spite several ba- 
sic differen c es in the respective models. Firstly, 
Goldsmith (l200ll) calculated the dust tempera- 
ture for a constant shielding that, in the absence 
of gas-dust coupling resulted in a temperature of 
Tdust ~6K. In our calculations the dust tempera- 
ture was solved self-consistently. The model was 
assumed to be illuminated by the full ISRF, and 
the radiation field inside the cloud was solved with 



radiative transfer calculations (|Juvelall2005l ). The 
difference in dust temperature has no effect on T gas 
at low ra(H 2 ). At n(H 2 ) = 10 5 cm Jrhi st. is close 
to the value used by I Goldsmith! ( 20011 ) and the 
T gas values are in very good agreement. Secondly, 
in our model the excitation temperature T ex varies 
as a function of radius so that the LVG assump- 
tion of a uniform medium is not valid. This affects 
the radiative connection between different parts of 
the model, is reflected in the excitation, and can 
thus affect the cooling rates. Thirdly, the photon 
escape probabilities in the LVG model and our 
Monte Carlo model are not identical even when 
the total optical depths are equal. In our calcula- 
tions, even at high optical depths, some photons 
can always escape in the line wings. These effects 
do not appear important for the estimated cen- 
tral temperature of the cores. There are also some 
differences. For example, in our calculations the 
gas-dust coupling has a smaller effect on T(j us t at 
high densities, apparently because our dust cool- 
i ng ra te is higher than that given by iGoldsmith 
(2001) Eq. 13. 



We used the homogeneous models to inves- 
tigate the effect of the grain size distribution. 
The gas-dust coupling becomes significant around 
rt(H 2 ) ~ 10 5 cm~ 3 depending, however, on the to- 
tal grain area. Equation. [T] is valid for a size dis- 
tribution dn/da ~ a~ 3,5 with a in the range 0.01- 
1.0 mm. If the low er limit is reduced to 10A (e.g. 
Li fc Drainell200lh . without modifying the gas-to- 



dust ratio, Kg.d increases by a factor of ^3. On 
the other hand, at the centre of dense cores the 
size of the large grains increases through grain 
coagulatio n while small grains may disappear en- 



tirely (e.g. IStepnik et al.l l2003l: lOrmel et al.ll200l 



Steinacker et al.ll2010l ). If the lower limit of grain 
sizes increases to 500A, the rate A g ^ is reduced 
by 60%. If the upper limit is further increased to 
2 /im, the effect is a factor of three. The value of 
A 9t d is similarly increased (decreased) by a factor 
of three if the powerlaw exponent of the size distri- 
bution is decreased (increased) by ~ 0.75, without 
modifying the size limits. 

Figure [2] illustrates the consequences for the gas 
temperature. The solid lines show the radial pro- 
files of Tgas and Tdust for the model n = 10 5 cm~ 3 
of Fig. [T] The increased photon escape proba- 
bility always decreases T gas at the cloud surface 
in spite of the increasing Td us t (Eq. [!}. The 
other curves correspond to a three times stronger 
and a three times weaker gas-dust coupling. The 
dash-dotted line is schematically the expected be- 
haviour where, for given density, the coupling be- 
comes weaker in the central part because of the 
grain growth. This is not to be taken so much as 
a model of an actual core as an illustration of the 
uncertainty of Tg as . 

3.2. Bonnor-Ebert spheres 

As more realistic models of dense cores, we 
examine critically stable 0.5 and 1.0 solar mass 
Bonnor-Ebert spheres. We will first investigate 
different factors that could affect their radial tem- 
perature profiles and then have a look at the pre- 
dicted line profiles. 

3.2.1. Temperature profiles 

We investigate first the one solar mass Bonnor- 
Ebert sphere. The basic model has a constant 
turbulent line width of cry = lkms -1 and no 
large scale gas motions. The central density rises 
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above 10 5 cm -3 but, because of the smaller cloud 
size, the column densities are lower than in the 
Sect. 13.11 models of similar density. The pho- 
ton escape probability increases outwards and the 
model predicts a significant decrease of T gas to- 
war ds_th^_cloud surface. This is in contrast with 
the I G alii et all (|2002t ) calculations that employed 
a parameterization of LVG results to estimate A ffas 
(see their Fig. 3; note also the difference in Td us t 
due to a different dust model) . The actual surface 
temperature will be sensitive to the abundance 
profiles and will depend on the amount of UV and 
cooling line radiation entering the cloud from the 
outside, both effects ignored in the present model. 

Figure [3] shows the quantitative effects result- 
ing from possible variations of the abundances, the 
velocity field, and A g ^. T he question of dep letion 
was already examined by iGoldsmithl (|200ll ) . To 
illustrate the effect in the context of our model, 
we decrease the abundance of all cooling species 
by a factor of ten in the cloud centre, within a ra- 
dius of 0.025 pc (Fig. [3^). Because many lines are 
already optically thick, A g is not expected to de- 
crease linearly with the abundances. In the model, 
T gas increases by ^2 degrees in the inner part, with 
only a small effect reflected in the outer cloud. 

In quiescent cores the line widths are some- 
times observed to be close to that determine d 
by thermal broadening (e.g. Hariu et al. 2008h . 



When the turbulent line width is reduced to ay = 
0.1 km s -1 within innermost 0.025 pc, the central 
temperature again increases by about two degrees 
(Fig. [3}d). The observed linewidths will remain 
much broader because of the thermal and opacity 
broadening and because of the emission from the 
outer cloud layers where the turbulent line width 
is still ay = 1.0 km s -1 . 

Although the Bonnor-Ebert model is static, we 
can introduce a large scale velocity field to check 
its importance on the escape of line emission. We 
add an infall velocity that is zero at the cloud sur- 
face and increases linearly to lkms -1 in the cen- 
tre. This is a very simplistic model of the velocity 
field but the m agnitude of the v elocity gradient 
is realistic (e.g. Zhou et al. 19931 ) and the model 
should capture the main effect on the radiative 
transfer. In the very central part of the model, 
the gas temperature by reduced little less than one 
degree (Fig. |3h). Both the large scale and small 
scale velocity fields affect A g through the line op- 



tical depths. 

The density of the one solar mass model ranges 
from 1.4 xlO 4 to 2xl0 5 cm 3 , a region of densities 
where the coupling between gas and dust becomes 
important. As discussed in Sect. 13.11 changes in 
the grain size distribution are reflected on the ef- 
ficiency of Kg : d so that it could be decreased by 
up to a factor of three in the cloud centre and en- 
hanced by up to a factor of three at the surface. 
It may be very improbable that a single source 
would exhibit the full range of variation. How- 
ever, Fig. [3ji shows this case where the efficiency 
of A 9i d jumps from one extreme to the other again 
at 0.025 pc radius. Because the dust temperature 
is above the gas temperature, the weaker coupling 
reduces T gas in the centre of the cloud. The effect 
is again of the order of one degree. 

Figure 0] shows the corresponding results for a 
half solar mass cloud where the central density is 
four times the value of the previous model. The 
main differences result from the enhanced gas-dust 
coupling that, together with lower dust tempera- 
ture, reduces the central gas temperature by ~ 1 K 
in the cases of low abundances and low velocity 
dispersion. On the other hand, the outer cloud 
is warmer by a similar amount, both because of 
the stronger A g> d but also because of the general 
density dependence already seen in Fig. Q] 

An external UV field can directly impact the 
cloud temperature at least at its surface. We 
examined another set of one solar mass Bonnor- 
Ebert spheres where the photoelectric heating was 
included and calculated as described in Sect. [2] 
The ISRF impinging on the cloud is attenuated 
by Ay=0, 1, or 2 m , corresponding to a shielding 
dust layer that is thought to exist around the ac- 
tual model cloud. This applies to the calculation 
of both Tdust and Tpe . The resulting temperature 
profiles inside the model are shown in Fig. [5j The 
effect of photoelectric heating becomes negligible 
once the cloud is shielded by Ay ~ 2 m . However, 
in an unshielded cloud the UV field has a small 
effect, ^0.4 K, even in the cloud centre. This is 
caused not by a direct photoelectric heating but 
by the change of the excitation in the outer cloud 
layers (cf. Sect. 0] Fig. [Tb). Again, the actual ef- 
fect will depend critically on the molecular abun- 
dances in the region heated by Tpe. 

As a final source of uncertainty we consider the 
rate of cosmic rays. The heating term T gaSyCr is 
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3x io-i7o-i 



based on the assumption of a rate £ = 3xiu~"' s 
but, in diffuse clouds, there are reports of rates 
that are up to two orders of magnitude higher (e 
McCall et all 120031: lLisztll2003tlShaw et al.l(200 



To reflect this uncertainty we calculated tempera- 
ture profiles with ( scaled by factors 1, 2, 5, and 
10. Figure [6] shows the resulting temperature pro- 
files, again for the one solar mass model. With 
the highest rate, ( = 3 x 10~ 16 s _1 , the central 
temperature has risen from the original ~7.5K to 
~17 K. For comparison, we show temperatures cal- 
culated with the A 90S parameterization given by 
I Goldsmith! ( 2001 ). The calculation is done shell 
by shell using the local density and the local dust 
temperature (r cr and A gaSi dust are as in our calcu- 
lations). These A gQS rates correspond to a model 
that has much higher column density (per velocity 
interval and for given density) than the Bonnor- 
Ebert spheres. Therefore, also the derived T gas 
values are higher. Furthermore, the parameteri- 
zation does not catch the increased photon escape 
probability at the cloud surface that, in our Monte 
Carlo calculations, results in the decrease of tem- 
perature in the outer part. 

3.2.2. Spectral lines 

Figures and HI show 13 CO and C 18 line pro- 
files that were calculated for the models of Figs. [3] 
andd The J =1-0, J =2-1, and J =3-2 spectra 
were calculated as observed towards the centre of 
the cloud with a beam with the FWHM equal to 
the half of the cloud radius. At 100 pc distance, 
this corresponds to ^52"and ^26"for the 1.0 M 
and the 0.5 M models, respectively. 

In the one solar mass model, the 13 CO(l- 
0) beam averaged optical depth is only r ^2.5 
through the cloud. In the 0.5 M Q the correspond- 
ing optical depth is ~7.5 meaning that there the 
13 CO transitions and C 18 lines originate par- 
tially in different regions with different kinetic 
temperatures. One example of this are the line 
profiles of the 0.5 M cloud with the infall veloc- 
ity. The 13 CO(2-l) line shows the expected in- 
fall profile while in the optically thinner J =1- 
line (r=3.9 vs r=6.5 for the second transition) 
the effect is weaker. The beam averaged C 18 0(1- 
0) lines remain symmetric. When observed with 
a pencil beam, the C 18 J=l-0 and J =2-1 lines 
would show slight asymmetry but with line pro- 
files with stronger emission on the red shifted side. 



This is caused by the T gas which decreases towards 
the cloud centre (see Fig. [4]:). 

We carried out LTE analysis of the 13 CO and 
C 18 lines to check how accurate those column 
density estim ates would be. We used the method 
described bv lMvers et al. I (ll983h . The ratio of the 
13 CO and C 18 lines is used to calculate the op- 
tical depth T\8 , and the excitation temperature of 
C 18 0, T e 18 , is solved from the radiative transfer 
equation (see iMvers et al.1 Il983l Eqs. 3-4). As- 



suming Gaussian line shapes, the column density 
of C 18 in the J = 1 state is 



N 



j=i 



3.86 x 10 14 r 18 J(T e 18 )Au 18 cm" 



(3) 



where J(T) = T Q /[exp(T /T) - 1} and T = 5.27 
K. The total C 18 column density is obtained by 
summing all levels, assuming they are populated 
according to T e 18 . Table [2] summarizes the results 
when the line parameters arc taken from Gaussian 
fits to the modelled spectral profiles. The errors 
of these estimates are less than 5% which shows 
that the temperature gradients have little impact 
on the column density estimates. 

4. The importance of non-local radiative 
couplings 

The LVG method is based on the assumption 
that the excitation is constant within the radia- 
tively coupled volume. This is not a good ap- 
proximation in dense cores where the gas veloc- 
ities are small and the radial gradients of T ex are 
large. To illustrate the potential problem fur- 
ther in a schematic way, we took from Fig. [T] the 
n = 10 5 cm~ 3 model and forced the inner part, 
r < 0.5 pc, to local thermodynamic equilibrium 
(LTE) at 5 K or 20 K. Fig. M shows the result- 
ing Tg as values in the outer part of the cloud, 
r > 0.5 pc. The assumed excitation of the inner 
cloud has a strong impact on the remaining cloud 
volume. The same effect was seen in Sect. 13. 2\ 
where the increased surface temperature caused 
by photoelectric heating was reflected all the way 
to the centre of the optically thick cloud. 

Our final example of a system with non-trivial 
radiative couplings consists of two Bonnor-Ebert 
spheres with properties identical to those in 
Sect. 13.21 The spheres are touching each other 
and the mutual shielding and exchange of ra- 
diative energy modifies the temperature distribu- 
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tions. Figure [TOl shows the T gas estimated without 
and with photoelectric heating. The calculations 
were carried out with the same continuum and 
line radiative transfer programs as in the case of 
spherical models but discretizing the cloud onto a 
128 3 cartesian grid. When photoelectric heating 
is not considered, the absorption of line radiation 
from the other core increases the temperatures 
between the cores. The effect on T gas is ~ 1 de- 
gree at the surface between the two cores. When 
UV heating is included, the mutual shielding be- 
comes important and T gas is reduced by up to 
3K between the cores. Compared to the gravi- 
tational attraction between the spheres, the force 
exerted by the resulting pressure asymmetry is 
of the order of one percent (assuming the change 
of AT ~ 3K affects a few percent of the surface 
area). However, if the spheres were to partially 
coalesce, the affected area and the temperature 
asymmetry would both increase making the effect 
potentially even dynamically important. 

5. Discussion 

We have modelled the gas and dust tempera- 
ture of dense clouds. The results emphasize the 
difference between the T gas distributions obtained 
with consistent radiative transfer calculation and 
those resulting from the blind application of LVG 
model results. If LVG calculations are used, also 
the increased photon escape probability near cloud 
surface must be taken into account. The effect can 
be several degrees and this will have consequences 
for chemical models and the interpretation of ob- 
servations. The photoelectric heating is capable of 
raising the surface temperature significantly but 
only if the cloud is shielded by less than Ay ~ l m 
of extinction. However, the indirect impact of this 
heating is felt well beyond the region directly pen- 
etrated by UV photons. The details of the radial 
temperature profile will depend on the abundances 
of the outer cloud layers and the external radiation 
field, both in terms of line radiation and the UV 
flux. Further studies coupling the chemistry and 
the modelling of the thermal balance are clearly 
needed. 

We examined the effects on T gas resulting from 
such spatial variation of molecular abundances, 
velocity field, and dust grain size distribution that 
are expected in dense cores. Each factor alone can 



change T gas by ~1 K or more. In the core, a strong 
depletion of molecules and the reduction of turbu- 
lent motions is capable of rising the t emperature 
by sev eral degrees. As pointed out by iGoldsmith 
([20011 ) . the effect of depletion will become less im- 
portant at higher densities when Ag^ dominates 
over line cooling. The same applies to any ef- 
fect resulting from the velocity field. However, 
the increase of grain sizes will significantly de- 
crease the coupling between gas and dust. When 
Idust > 7g as and the density is close to 10 5 cm" 3 , 
this can compensate some of the temperature in- 
crease predicted for the inner core. There is no 
observational evidence of a temperature increase 
at the centre of quiescent cores but also this possi- 
bility should be considered when interpreting ob- 
servations. In more opaque clouds (especially in 
conjunction with dust coagulation) Td us t will be 
reduced below T gas and, at high enough densi- 
ties, will eventually force gas temperature down 
at the centre of starless cores. An increase in 
the grain sizes can shift this transition to densi- 
ties higher than usually assumed. However, be- 
cause of_the_l_ong time scale of dust coagulation 
(e.g. lOrmel etaA lioooh . the effect is often likely 
to be smaller than in our model. 

The value of T gas is particularly uncertain at 
the cloud surface where the photoelectric heating 
and the dissociation of CO produce strong tem- 
perature gradients. However, for the dense clouds 
the largest source of uncertainty still appears to be 
the rate of cosmic ray heating. In the case of the 
one solar mass Bonnor-Ebert sphere, a factor of 
five increase in £ would raise the central tempera- 
ture by five degrees to T gas ^12.5 K. This may al- 
ready be excluded by direct observational evidence 
of much lower gas temperatures in d ense clouds 
(|Crapsi et alj|2007t Irlariu et al.ll2008l) . Neverthe- 
less, the theoretical prediction of the temperatures 
- and temperature profiles - of dense cores still 
contains significant uncertainty. 

The observed temperature gradients will not 
strongly modify the radial density distribution of 
cores nor significantly affect the cor e stability (e.g. 
Hariu et al.fl2008t iGalli et al.l l2002h . However, in 



the critically stable Bonnor-Ebert models of one 
solar mass, the difference between isothermal tem- 
peratures of 8K and 10 K corresponds to a fac- 
tor of two increase in the central density. This 
shows that even small temperature changes are 
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important in theoretical studies. In Figs. [3] and 
2] we examined separately the effect of various 
parameters on T gas . In more dense cores, with 
Tdust < the effects of the x, cr v , and A g ^ d pa- 
rameters could accumulate, making the tempera- 
ture gradients more pronounced. If gas is not cou- 
pled to dust, T gas could in the central parts remain 
several degrees above the temperature of the outer 
cloud. However, in the centre of the 0.5 Mq model 
cloud the gas temperature was already largely de- 
termined by the gas-dust coupling. In that case 
the question of the dust properties becomes im- 
portant because different dust models can, in the 
cloud centre, lead to Tdust values differi ng by more 
than one degree (e.g. our Td us t vs. the (I G alii et al 



2002h models). The case of the two Bonnor-Ebert 



spheres in Sect. [4] suggests that even small temper- 
ature anisotropics may sometimes play a role in 
the long term evolution of clouds. Further studies 
are also required to find out how the early evolu- 
tion of spherical cores is modified relative to the 
isothermal case. 

However, the main importance of small T gas 
variations may come via chemistry. The gas tem- 
perature directly affects chemical reaction rates 
and, in particular, the depletion onto dust grains. 
Therefore, the precise value of the gas tempera- 
ture is relevant for interpretation of both line and 
continuum data. Because the collision rates are 
only proportional to y/T gas , the direct effect on 
the time scales of depletion and grain mantle ac- 
cumulation is small, less than 50% for the kind 
of temperature variations observed in our models. 
However, the effect on steady state abu ndances 



is more noticeable. lAikawa et al.l ( 20051 ) studied 



chemical evolution in collapsing clouds with initial 
conditions close to critical Bonnor-Ebert spheres. 
If the evolution was slow enough, significant de- 
pletion was observed in the central parts of the 
model clouds. In particular, Aikawa et al.l ( 2005 ) 
included a comparison of identical models (central 
density 3 x 10 6 cm~ 3 ) with kinetic temperatures 
of 10 K, 12 K, and 15 K. The ice composition was 
found to be very sensitive to the temperature and 
this was reflected in the gas phase abundances. A 
difference of two degrees could modify some abun- 
dances by a factor of two and, in the centre where 
the depletion becomes significant, by an order of 
magnitude or even more (e.g. NH 3 and N 2 H + , see 



Aikawa et alj 1200a , Figs. 2 and 6). It is conceiv- 



able that in some cases the depletion will be reg- 
ulated by the temperature rise that results from 
the decreasing line cooling. 

When observations are analyzed, there is no 
guarantee that different lines (e.g., different iso- 
topomers or different transitions of the same 
molecule) would originate in identical gas vol- 
ume. We already noted that, depending on the 
opacity of the lines, the spatial resolution, and 
the radial T gas profile, it is possible to observe 
both blueshifted and redshifted spectra towards 
a collapsing cloud. In the same fashion, the ki- 
netic temperature measured, e.g. with NH3, may 
not be representative for other lines. This could 
lead to errors that are propagated to the derived 
column densities. The evaluation of these un- 
certainties requires simultaneous modelling of the 
thermal balance, radiative transfer and, in par- 
ticular, of the chemistry. Such a full study is 
beyond the scope of the present paper. However, 
we did carry out LTE analysis of the 13 CO and 
C 18 lines calculated for the Bonnor-Ebert mod- 
els. The temperature gradients did not affect the 
column density estimates by more than 5% and, at 
least in this case, the LTE analysis would produce 
accurate estimates for the total 13 CO and C 18 
column density. Thus, the temperature variations 
examined in this paper mainly affect our expecta- 
tions of the radial abundance profiles. These are 
important considerations when line data are used 
to estimate the central density or temperature of a 
dense core. The direct implications on the stabil- 
ity or dynamic evolution of the cores are probably 
of secondary importance. 

MJ and NY acknowledge the financial support 
by the Academy of Finland Grant 127015. 
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Table 1 

The list of modified parameters affecting T e 



Parameter Default value Range of values" 

A s ,d Eq.Q] 0.3-3.0 x A° d 

abundances x Goldsmith (2001) Table 1 0.1-1.0 xx° 

infall velocity Okms -1 0.0-l.Okms" 1 

turbulent linewidth a v l.Okms - 1 0.1— l.Okms -1 

rate of cosmic rays C 3xl0 _17 s _1 1.0-10.0 xC° 



a Thc default parameter value is indicated with an upper index 
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4 5 

MH 2 )/cm 3 ) 



Fig. 1. — A comparison of the T gas (squares) and 
Tdust ( circles) in o u r calc ulations (large symbols) 
and in iGoldsmithl (|200ll ) (small symbols). The 
open symbols denote the results without the gas- 
dust coupling. 
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Fig. 2. — The radial profiles of Tdust (lower solid 
line) and T gas (upper solid line) for a homogeneous 
cloud with a density of n = 10 5 cm~ 3 . The other 
lines correspond to three times stronger or three 
times weaker coupling or a combination of the two 
where the coupling is weaker in the centre. 
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R (pc) 
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R (pc) 



Fig. 3. — Tg as in a one solar mass Bonnor-Ebert 
sphere. The dashed line shows the default solution 
with constant fractional abundances x, a turbu- 
lent line width of cry = 1 km s~ 1 , and no large scale 
velocity gradients. Part of the dust temperature 
profile is shown in the first frame (Td us t >11K). 
Each frame shows the effect of a single modifica- 
tion: a lower fractional abundance for r < 0.025 pc 
(frame a), lower turbulence try = 0.1 kms -1 for 
r < 0.025 pc (frame b), an infall velocity (frame 
c), and a lower value of A g ^ (frame d). See text 
for details. 
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Fig. 4, 
(lines as in Fig. |3]) 



Tg as in 0.5 Mp, Bonnor-Ebert models 
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(f.00 0.01 0.02 0.03 0.04 0.05 0.06 

R (pc) 

Fig. 5. — The radial profiles of T gas (solid lines) 
and T(j us t (dashed lines) for the one solar mass 
Bonnor-Ebcrt model, taking into account the pho- 
toelectric heating. Reading from the top, the lines 
correspond to clouds shielded by an Ay equal to 
0, 1, or 2 magnitudes. 
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b. x 


_ C. f rad 


d. A(gas—dust) 



-3-2-10 1 2 3-3-2-10 1 2 3 



V{hn/s) V(km/s) 

Fig. 7. — The 13 CO spectra (solid lines) and the 
C 18 spectra (dashed lines) for selected models 
from Fig. |3] The curves correspond to transitions 
,7 = 1 — (thick lines), J = 2 — 1 (medium line), 
and J = 3 — 2 (thin lines). The C ls O spectra have 
been scaled by a factor of two. 
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Fig. 6.— The radial profiles of T gas in 1.0 Af© 
Bonnor-Ebert models where the cosmic ray heat- 
ing rg as . cr is scaled by factors 1.0, 2.0, 5.0, and 
10.0 (the curves reading from the bottom). The 
dashed line s are te mperatures calculated using the 
Goldsmith! (|200lh parameterization of A gas that 
applies to LVG models with higher column den- 
sity. 
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Fig. 8.— The 13 CO and C 18 spectra for the 
0.5 M Q models of Fig.0] The lines are as in Fig. [7] 
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Table 2: LTE column density estimates derived for 
the model spectra of Figs. [7]and|Hl 



Model version 


1.0 Mq 


model 


0.5 M Q 


model 




^V(C 18 0)? ruo 


jV(C 18 0) LTE 


^(C 18 0)f ruc 


iV(C 18 0) LTE 




(cm 2 ) 


(cm -2 ) 


(cm 2 ) 


(cm -2 ) 


Default 


1.79xl0 15 


1.88xl0 15 


3.57xlO ib 


3.75xlO ib 


Modified \ 


6.72xl0 14 


6.97xl0 14 


1.34xl0 15 


1.38xl0 15 


Modified w ra d 


1.79xl0 15 


1.88xl0 15 


3.57xl0 15 


3.74xl0 15 


Modified A g ,d 


1.79xl0 15 


1.89xl0 15 


3.57xl0 15 


3.74xl0 15 



a The true, beam averaged column density of the model 
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Fig. 9. 


- For the n = 10 5 cm~ 3 


model, the radia. 



temperature profile in the outer part of the cloud. 
The solid curve is the solution from Fig. [TJ The 
dashed line and the dotted line indicate the tem- 
perature when the inner part of the core is forced 
to LTE with a temperature of 20 K or 5 K, respec- 
tively. 
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Fig. 10. — Gas temperature in a system consist- 
ing of two one solar mass Bonnor-Ebert spheres 
that are touching each other. The images on the 
top show T gas in a cross section through the cen- 
tre of the lower sphere, the symmetry axis be- 
ing vertical. The left and right hand sides corre- 
spond, respectively, to the situation without and 
with photoelectric heating. The lower frame shows 
the surface temperature, averaged over 90-95% of 
the core radius as the function of the angle 9. 
The angle is measured from the symmetry axis, 
as seen from the centre of a core, and the value of 
6 = 0° corresponds to the direction towards the 
other sphere. The dashed line and the right hand 
axis show T gas with photoelectric heating included. 



13 



